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Abstract —In this paper, a new mathematical framework to 
the analysis of millimeter wave cellular networks is intro¬ 
duced. Its peculiarity lies in considering realistic path-loss and 
blockage models, which are derived from experimental data 
recently reported in the literature. The path-loss model accounts 
for different distributions for line-of-sight and non-line-of-sight 
propagation conditions and the blockage model includes an 
outage state that provides a better representation of the outage 
possibilities of millimeter wave communications. By modeling the 
locations of the base stations as points of a Poisson point process 
and by relying upon a noise-limited approximation for typical 
millimeter wave network deployments, exact Integral expressions 
for computing the coverage probability and the average rate 
are obtained. With the aid of Monte Carlo simulations, the 
noise-limited approximation is shown to be sufficiently accurate 
for typical network densities. Furthermore, it is shown that 
sufficiently dense millimeter wave cellular networks are capable 
of outperforming micro wave cellular networks, both in terms of 
coverage probability and average rate. 

I. Introduction 

In spite of common belief, recently conducted channel 
measurements have shown that millimeter wave (mmWave) 
frequencies may be suitable for cellular communications, 
provided that the cell radius is of the order of 100-200 
meters [1]. Based on these measurements, the authors of [2] 
have recently studied system-level performance of mmWave 
cellular networks and have compared them against conven¬ 
tional micro wave (/rWave) cellular networks. The results 
have highlighted that mmWave cellular communications may 
outperform /rWave cellular communications, by assuming 
similar cellular network densities, provided that a sufficient 
beamforming gain is guaranteed between Base Stations (BSs) 
and Mobile Terminals (MTs). These preliminary but encour¬ 
aging results have motivated several researchers to investigate 
potential and challenges of mmWave cellular communications 
for wireless access, in light of the large and unused bandwidth 
that is available at these frequencies. 

System-level performance evaluation of cellular networks is 
widely recognized to be a mathematically intractable problem 
[3], due to the lack of tractable approaches for modeling 
the other-cell interference. Only recently, a new mathematical 
methodology has gained prominence due to its analytical 
tractability and reasonable accuracy compared to currently 
deployed cellular networks. This emerging approach exploits 
results from stochastic geometry and relies upon modeling the 
locations of the BSs as points of a point process [3]. Usually, 


the Poisson Point Process (PPP) is used due to its mathematical 
tractability [4]. Recent results for cellular networks modeling 
based on stochastic geometry are available in [5]-[9], to which 
the reader is referred for further information. 

Motivated by the mathematical flexibility of the PPP-based 
abstraction modeling, researchers have recently tried to study 
the performance of mmWave cellular communications with 
the aid of stochastic geometry, by developing new frame¬ 
works specifically tailored to account for the peculiarities of 
mmWave propagation [10], [11]. In fact, the mathematical 
frameworks currently available for modeling /iWave cellular 
networks are not directly applicable to mmWave cellular 
networks. The main reasons are related to the need of incor¬ 
porating realistic path-loss and blockage models, which are 
significantly different from /rWave communications. Eor ex¬ 
ample, [1] and [2] have pointed out that Line-Of-Sight (LOS) 
and Non-Line-Of-Sight (NLOS) links need to be appropriately 
modeled and may have different distributions, due to the more 
prominent impact of spatial blockages at mmWave frequencies 
compared to /rWave frequencies. Eurthermore, in mmWave 
communications a new outage state may be present in addition 
to LOS and NLOS states, which better reflects blockage effects 
at high frequencies and accounts for the fact the a link may 
be too weak to be established. 

Recently reported results on stochastic geometry modeling 
of mmWave cellular communications [10], [11] take these as¬ 
pects into account only in part. The approach proposed in [10] 
is mathematically tractable for dense network deployments 
and relies upon an equivalent LOS ball approximation. The 
approach proposed in [11] uses a similar LOS ball approxima¬ 
tion, but is applicable to medium/sparse network deployments. 
The outage state observed in [2] is not explicitly taken into 
account either in [10] or in [11]. Against this background, 
in the present paper a new methodology to the stochastic 
geometry modeling of mmWave cellular communications is 
proposed, which explicitly incorporates the empirical path-loss 
and blockage models recently reported in [2]. 

The paper is organized as follows. In Section [III system 
model and main modeling assumptions are introduced. In Sec¬ 
tion m the mathematical frameworks for computing coverage 
probability and average rate are described. In Section IIVI the 
frameworks are validated via numerical simulations and the 
performance of mmWave and /rWave cellular networks are 
compared. Einally, Section IV] concludes this paper. 


II. System Model 

A. PPP-Based Abstraction Modeling 

A bi-dimensional downlink cellular network deployment is 
considered, where a probe MT is located, without loss of 
generality thanks to the Slivnyak theorem [12, vol. 1, Th. 
1.4.5], at the origin and the BSs are modeled as points of 
a homogeneous PPP, denoted by 'P, of density A. The MT is 
assumed to be served by the BS providing the smallest path- 
loss to it. The path-loss model is introduced in Section III-DI 
The serving BS is denoted by Similar to [3, Sec. VI], 

full-frequency reuse is considered. For notational simplicity, 
the set of interfering BSs is denoted by = 'P\BS*'°^ 

The distance from the generic BS to the MT is denoted by r. 


B. Directional Beamforming Modeling 

Thanks to the small wavelength, mmWave cellular networks 
are capable of exploiting directional beamforming for com¬ 
pensating for the increased path-loss at mmWave frequencies 
and for overcoming the additional noise due to the large 
transmission bandwidth. As a desirable bonus, directional 
beamforming provides interference isolation, which reduces 
the impact of the other-cell interference. Thus, antenna arrays 
are assumed at both the BSs and the MT for performing 
directional beamforming. For mathematical tractability and 
similar to [10] and [11], the actual antenna array patterns 
are approximated by considering a sectored antenna model. 
In particular, the antenna gain of the generic BS, Gbs (•)’ 
of the MT, Gmt (•)’ ^e formulated as follows: 


Gq(0) 


^(max) 

g(“”) if |0|>u;q 


( 1 ) 


where q € {BS,MT}, 0 € [0,27r) denotes the angle off the 
boresight direction, ujq denotes the beamwidth of the main 
lobe, is the array gain of the main lobe and G,™™^ is 

the antenna gain of the side lobe. 

Both the MT and its serving BS, BS^°\ are assumed to 
estimate the angles of arrival and to adjust their antenna steer¬ 
ing orientations accordingly. On the intended links, therefore, 
the maximum directivity gain can be exploited. Thus, the 
directivity gain of the intended link is G*^°^ = . 

The beams of all non-intended links are assumed to be ran¬ 
domly oriented with respect to each other and to be uniformly 
distributed in [0, 27r). Accordingly, the directivity gains of the 
interfering links, G^®^ for i G 'pGo)^ are randomly distributed. 
Based on ([T]i, their Probability Density Function (PDF) can be 
formulated as fcn) {g) = (f a) (g) f a) ) (g), where: 

\ ^BS ^MT/ 

‘s’=S'* (’ - °»"“’) + - S) '* (»- 

( 2 ) 

for q G {BS, MT}, S (•) is the Kronecker’s delta function and 
g) denotes the convolution operator. 


C. Link State Modeling 

Let an arbitrary link of length r, i.e., the distance from 
the generic BS to the MT is equal to r. Motivated by recent 
experimental findings on mmWave channel modeling [2, Sec. 
III-D], a three-state statistical model for each link is assumed, 
according to which each link can be in a LOS, NLOS or in 


an outage (OUT) state. A LOS state occurs if there is no 
blockage between BS and MT. A NLOS state, on the other 
hand, occurs if the BS-to-MT link is blocked. An outage state 
occurs if the path-loss between BS and MT is so high that 
no link between them can be established. In this case, the 
path-loss of the link is assumed to be infinite. In practice, 
outages occur implicitly when the path-loss in either a LOS 
or a NLOS state are sufficiently large. In [2, Fig. 7], with the 
aid of experimental results, it is proved that adding an outage 
state, which usually is not observed for transmission at /rWave 
frequencies, provides a more accurate statistical description of 
the inherent coverage possibilities at mmWave frequencies. 

In particular, the probabilities of occurrence plos (•)’ 
Pnlos (•)’ POUT (•) of LOS, NLOS and outage states, respec¬ 
tively, as a function of the distance r can be formulated, based 
on experimental results, as follows [2, Eq. 8]: 

POUT (t) = max (0,1 - 7ouTe“'^°'"'^’'} 

PLOS (t) = (1 - POUT (t)) 7LOse“'^’^°®'’ (3) 

Pnlos (t) = (1 - POUT (t)) (1 - 7LOse"'^^°®’') 

where (^los, 7los) and ((5 out, 70 ut) are parameters that de¬ 
pend on the propagation scenario and on the carrier frequency 
being considered. Examples are available in [2, Table I]. 

Under the assumptions that the BSs are modeled as points of 
a homogeneous PPP and that the probabilities of the BS-to-MT 
links being in a LOS, NLOS or outage state are independent, 
4/ can be partitioned into three (one for each link state) 
independent and non-homogeneous PPPs, i.e., 'Flos^ 'Fnlos 
and T'out? such that 'F = 'Flos O ^^nlos O 'Fqut- This 
originates from the thinning property of the PPPs [12]. Prom 
©, the densities of the PPPs 'Flos^ 'Fnlos T'out are 
equal to Alos (t) = ApLOS (t), Anlos (t) = Apnlos (t) and 
Aout (t) = ApouT (t), respectively. 

D. Path-Loss Modeling 

Based on the channel measurements reported in [2], the 
path-loss model for LOS and NLOS links is assumed to be: 

^Los (t) = , Wos (t) = (4) 

where r denotes the generic BS-to-MT distance, klos and 
knlos can be interpreted as the path-loss of LOS and NLOS 
links at a distance of 1 meter, respectively, /3 los and /3 nlos 
denote the power path-loss exponents of LOS and NLOS links, 
respectively. 

The path-loss model in (|4|i is general enough for modeling 
several practical propagation conditions. Por example, it can be 
linked to the widespread used (a, /3) model [1], [2], by setting 
klos = 10“''°s/(io/3los) and knlos = 10“n''°s/(io/3nlos)^ 
where olos and aNUOS are defined in [2, Table I]. As 
mentioned in Section ITl-Cl the path-loss of the links that are in 
an outage state is assumed to be infinite, i.e., Zqut (t) = oo. 

E. Fading Modeling 

In addition to the distance-dependent path-loss model of 
Section III-Dl each link is assumed to be subject to a complex 
randomly distributed channel gain, which, for a generic BS- 
to-MT link, is denoted by h. According to [2], the power gain 
|/i|^ is assumed to follow a Log-Normal distribution with mean 
(in dB) equal to and standard deviation (in dB) equal 








Re-^) n{x- Zs) 
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to Thus, \hf' takes into account large-scale shadowing. 

In general, and for LOS and NLOS links are 

different [2]. In what follows, they are denoted by and 

where s = {LOS, NLOS} denotes the link state. 

As mentioned in Section ITl-CI for mathematical tractability, 
the shadowing correlations between links are ignored. Thus, 
the fading power gains of LOS and NLOS links are assumed 
to be independent and identically distributed. As recently 
remarked and verified with the aid of simulations in [10], this 
assumption usually causes a minor loss of accuracy in the 
evaluation of the statistics of the Signal-to-Interference-plus- 
Noise-Ratio (SINR). Due to space limitations and for ease of 
description, fast-fading is neglected in the present paper, but 
it may be readily incorporated. 

F. Cell Association Criterion 

The MT is assumed to be served by the BS providing the 
smallest path-loss to it. Let LL*og, -^out 

smallest path-loss of LOS, NLOS and OUT links, respectively. 
They can be formulated as; 

\+oo if = 0 (5) 

= min jZouT | = +00 

where s = {LOS, NLOS}, denotes the distance of a 
generic BS to the MT, and 0 denotes an empty set. Hence, 
the path-loss of the serving BS, BS^°\ can be formulated as 
r ( 0 ) _ / r ( 0 ) r (0) r (0) \ 

U — mm |ULQg, Uj^LOS’-'^OUT J- 

III. Mathematical Modeling of Coverage 
Probability and Average Rate 
Based on the system model of Section [III the SINR of the 
cellular network under analysis can be formulated as follows: 
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where P is the transmit power, | | ^ is the power gain of 


the serving BS, BS^°\ ajf is the noise power and Jagg is 
the aggregate other-cell interference, i.e., the total interference 
generated by the BSs in v[/(\o) in particular, cr^ is defined 
as a% = io'^«(dBm)/io^ ^^(dBm) = -174 + 


lOlogj^Q (BW)+J{iB, BW is the transmission bandwidth and 
Jtjb is the noise figure in dB. By using a notation similar to 
the desired link, the aggregate other-cell interference is defined 
as 4gg = Eietcvo) 1^^*^ r/■ The equality in (a) 

takes into account that the distribution of LOS and NLOS links 
is assumed to be different in the present paper. 

From coverage probability (Pcov) and average rate (R) 
can be formulated as follows [8]: 

p(cov) ^ Pj. {sink > T} (7) 

R = Esinr {BWlog2 (1 + SINR)} 

_ BW /■+°° Pcov (f) ,, (8) 

ln{2)Jo t+1 

where T is a reliability threshold and E{-} denotes the 
expectation operator. 

In what follows, a new mathematical expression for the 
coverage probability, (•), is provided. The average rate, 

R, can be obtained directly from ®. The analytical for¬ 
mulation is based on the noise-limited approximation of 
mmWave cellular communications, i.e., SINR « SNR = 
1^ (tT^L(°^) which has been observed in recent 
studies, both with the aid of numerical simulations and field 
measurements [2], [11]. In Section ITVl the validity and the 
accuracy of the noise-limited approximation are substantiated 
with the aid of Monte Carlo simulations, which account for 
the other-cell interference as well. 

Before presenting the main result summarized in Proposi¬ 
tion Q] we introduce two lemmas that are useful for its proof. 

Lemma 1: Let L = {Llos, Tnlos, Tout}, where Lg = 
{Zg (r^"^) ,n G T's} for s G {LOS, NLOS, OUT} are trans¬ 
formations of the path-loss of LOS, NLOS and OUT BSs, 
respectively. Let the link state and the path-loss models in Q 
and (IH), respectively. L is a PPP with intensity as follows: 

A ([0,x)) = Alos ([0,t:)) + Anlos ([0,a:)) (9) 

where ALos([0,a:)) = To(x;s = LOS), Anlos ([0, a:)) = 
Ti (x; s = NLOS)-To (x; s = NLOS), To (•; •) and Ti (•; •) 
are defined in Coll, R (■) is the Heaviside function, 
Rix) = 1 - n (x), /Cl _ = 27rA7LOS^Los, IC 2 = 

27rA7Los7ouT (<5 los + i^out) , R = ^los<5out (tout), 
W = (i5los + ^out) <5out (tout), Qs = <5loskJ^, 
Ts = SoUTK~^, Vs = (i5lOS + i^OUt) Zs = 

In (tout))^* for s = {LOS, NLOS}. 

Proof: See the Appendix. □ 

Corollary 1: Let i5out = 0, 70 UT = 1, i.e., 
POUT (r) = 0 in ©. A(-) in ® holds with 
To(a;;s) = ICi (l - 
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(a;;s) = (e"^ + We-^ - - Z,) 
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+ IC 2 {V^/Ps) {x - Zs) + /Cl (Q^Ps) {x - ZP) 

{x\ s) = -■nXnp'^x^/^^d {x - Zs) + ttA (^qut (tout))^ <5 {x - Zs) 

+ 27rA(5ouT70UT (tqot + 7oi1jt (tout) - - TsX ^/) <5 (a: - Zs) 

+ 2tt\kP^PP^x^/P‘-^H {x - Zs) + 2TT\5pl^-iovTT^PP^x^/^^-^e-^‘^"''‘n {x - Zs) 


Ti (x; s) = ttXkJ^x^/1^% s = {LOS, NLOS}. 

Proof: It follows directly from ( fTOl i. since Zs ^ +00 for 
s = {LOS, NLOS}. □ 

Lemma 2: Let L^°'> = min |L[°(^g, L^°^og, = 

min {Lj be the smallest element of the PPP L introduced 
in Lemma\I\ Its Cumulative Distribution Function (CDF), i.e., 
^L(o) i.^) = P'' < a;}, can be formulated as follows: 

Pl(o) (a:) = 1 - exp(-A([0,x))) (11) 

where A (•) is defined in ®. 

Proof: It follows by applying the void probability theorem 
of PPPs [13, Corollary 6 ], □ 

A. Coverage Probability 

Proposition 1: Let the SINR in (| 6 ]) and the coverage prob¬ 
ability in ©. If a% > Jagg, i-e., SINR r; SNR = 

then the following holds: 

P(“V) (T) = p[{°"g) (T) + P^^o^g (T) (12) 


p(cov) /rrix _ 1 / 

^ ^ ^ ^Jo 




sf2o 


As ([ 0 ,x))dx 

(13) 


where s = {LOS, NLOS}, erfc(-) is the complementary 
eiTor function, As([0,x)) = aI^^ ([0, x)) exp (—A ([0, x))), 
A(-) is defined in (|9]l, ALos([0ia:)) = Tg^^(x;LOS), 
where Tg^^(-;-) is the first derivative of To(-;') in (fTOl l. 
^nlos([ 0>2:)) = (x;NLOS) - (x;NLOS), where 
(•;•) is the first derivative of Ti (•;•) in (ITOl l. 7 !°! = 
ja\q, fis = (10)/10 and as = (10)/10. 

Tg^^ (•; •) and (•; •) are available in ([l4l i. 

Proof: See the Appendix. □ 


IV. Numerical and Simulation Results 

In this section, we illustrate some numerical examples 
for validating the accuracy of the proposed mathematical 
frameworks and for comparing mmWave and /rWave cellular 
networks. The results are compared against Monte Carlo 
simulations, where some modeling assumptions considered 
for developing the mathematical frameworks are not enforced 
in the system simulator. Notably, coverage probability and 
average rate in Section [III] are computed under the assumption 
of noise-limited cellular networks. This assumption is not 
retained in the system simulator, in order to show to which 
extent the noise-limited assumption holds for mmWave cellular 
communications. Monte Carlo simulation results are obtained 


by using the system simulator described in [5]-[8], to which 
the reader is referred for further information. 

Unless otherwise stated, the following setup is considered 
for obtaining the numerical examples, which agrees with 
previous studies in this field [2], [10], [11]. In particular, the 
channel and blockage models are taken from [2]. In addition: 

• Two mmWave cellular networks are studied, which oper¬ 
ate at a carrier frequency, F^, equal to = 28 GHz and 
Fc = 73 GHz. The transmission bandwidth is BW = 2 
GHz. The noise figure is A/ub = 10. The transmit power 
is P = 30 dBm. 

• The path-loss model is as follows [2, Table I]: aLOS = 
61.4 dB, /3los = 2 and onlos = 72 dB, /3nlos = 2.92 
if Fc = 28 GHz and cilos = 69.8 dB, /3los = 2 and 
aNLOS = 82.7 dB, /3nlos = 2.69 if F, = 73 GHz. 

■ The shadowing model is as follows [2, Table I]: crj^g^ = 

5.8, = 8.7 if F^ = 28 GHz and cr[‘Jfg^ = 5.8, 

®’ot!os “ 8.7 if Fc = 73 GHz. On the other hand, 

is assumed to be equal to zero for both LOS and NLOS 
scenarios. 

■ The blockage model is as follows [2, Table I]: i5los = 
1/67.1, 7los = 1 and (5out = 5.2, 70 UT = exp (1/30), 
for both Fc = 28 GHz and Fc = 73 GHz scenarios. 

• The directional beamforming model is as follows [10]: 

^BS “ ^MT “ ^BS “ ^MT ~ 

and wbs = = 30 degrees. 

• Similar to [10], the density of BSs, A, is represented as a 
function of the average cell radius, i.e., Rc = ^Ji/JttX). 

• As for the pWave cellular network, a setup similar to [2] 

is considered. In particular, we set Fc = 2.5 GHz, BW = 
40 MHz, = 0 dB and wmt = 360 

degrees. The channel model is chosen as in [2, Eq. (11)], 
i.e., I = 22.7-f 36.7logio (0 +261ogig (2.5). All 

the channels are assumed to be in a NLOS state, with a 
shadowing standard deviation equal to unlos = 4. No 
outage state is assumed, i.e., pout (t) = 0. The rest of 
the paraments is the same as for the mmWave cellular 
network setup. 

Some selected numerical results are illustrated in Figs. [T]- 
|4] From Figs. [T] and |2] we observe that the proposed noise- 
limited approximation is quite accurate for practical densities 
of BSs. If Rc > 100 meters for the considered setup, in 
particular, we observe that mmWave cellular networks can be 
assumed to be noise-limited. If the density of BSs increases, 
on the other hand, this approximation may no longer hold. 
The performance gap compared to Monte Carlo simulations is, 
however, tolerable and this shows that, in any case, mmWave 
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Fig. 1. Coverage probability of a mmWave cellular network operating at 
Fc = 28GHz. (a)pouT(')inO) is used, (b) pouxC^”) = 0. Solid lines show 
the mathematical framework and markers show Monte Carlo simulations. 
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Fig. 2. Coverage probability of a mmWave cellular network operating at 
Fc = 73 GHz. (a) Pout(') in ^3) is used, (b) pouT(f) = 0. Solid lines show 
the mathematical framework and markers show Monte Carlo simulations. 


cellular networks are likely not to be interference-limited. 
This finding is in agreement with recent published papers that 
considered a simplified blockage model [11]. The figures also 
show that, in general, the presence of an outage state reduces 
the coverage probability. This is noticeable, however, only for 
small values of the reliability threshold T. 

In Fig. [3 mmWave and /rWave cellular networks are com¬ 
pared. This figure shows that mmWave communications have 
the potential of outperforming /rWave communications, pro¬ 
vided that the network density is sufficiently high. Otherwise, 
/rWave communications are still to be preferred, especially 
for small values of the reliability threshold T. As expected, 
mmWave transmission at Fc = 28 GHz slightly outperforms 
its counterpart at Fc = 73 GHz due to a smaller path-loss. 

Finally, Fig. |4] shows the average rate of mmWave and 
/rWave cellular networks. It confirms that the proposed frame¬ 
work is accurate enough for practical network densities. Also, 
it highlights that mmWave communications are capable of 
significantly enhancing the average rate. This is mainly due to 
the larger transmission bandwidth, which is 50 times larger. 


Fig. 3. Coverage probability of mmWave and /iWave cellular networks 
operating at Fc = 28 GHz (mmWave), Fc = 73 GHz (mmWave) and Fc = 
2.5 GHz (pWave). For mmWave networks, Pout(') in G) is used. For the 
pWave cellular network, pouT(r) = 0. (a) Rc = 50 m. (b) Rc = 100 m. 
(c) Rc = 150 m. (d) Rc = 200 m. 


(a) (b) 



Fig. 4. Average rate of a mmWave cellular network operating at Fc = 
28 GHz and Fc = 73 GHz. (a) The normalized rate R/BW is shown. 
Solid lines show the mathematical framework and markers show Monte Carlo 
simulations, (b) Ratio of the average rates of two mmWave networks operating 
at Fc = 28 GHz and Fc = 73 GHz and of a /rWave network operating at 
Fc = 2.5 GHz. 


in the considered setup, for mmWave communications. The 
figure shows, however, that the gain can be larger than the 
ratio of the bandwidths, especially for medium/dense cellular 
network deployments. 

V. Conclusion 

In the present paper, a new analytical framework for com¬ 
puting coverage probability and average rate of mmWave cel¬ 
lular networks has been proposed. Its novelty lies in taking into 
account realistic channel and blockage models for mmWave 
propagation, which are based on empirical data recently re¬ 
ported in the literature. The proposed methodology relies upon 
the noise-limited assumption for mmWave communications, 
which is shown to be sufficiently accurate for typical densities 
of BSs. The numerical examples have shown that sufficiently 
dense mmWave cellular networks have the inherent capability 
of outperforming /rWave cellular networks. 
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Current research activities are concerned with the exten¬ 
sion of the proposed mathematical framework to ultra-dense 
mm Wave cellular networks, where interference can no longer 
be neglected, as well as to compare different blockage models 
and cell association criteria. An extended version of the present 
paper is available in [14]. 

Appendix - Proofs 


where erf (•) is the error function, (a) follows from the CDF of 
a Log-Normal random variable, and (b) and (c) follow from 
Lemma |7] and Lemma \2\ by using arguments similar to the 
computation of the intensity of L and of the CDF of L^^\ 


The proof follows by explicitly writing the expectation with 
respect to L^qq and L^los terms of their PDFs, which 
can be formulated, similar to (b) and (c), as f^(o) (0 = 

dPr ([0,^))exp(-As ([0,^))), since 

Pr |Li°^ < C} = exp (-As ([0,0)). 
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A. Proof of Lemma [7] 

The proof follows by using a methodology similar to [13, 
Sec. II-A]. In particular, by invoking the displacement theorem 
of PPPs [12, Th. 1.10], the process of the propagation losses 
L = can be interpreted as a transformation 

of which is still a PPP on M’*’. From Section [Tl-CI we know 
that ^ = ^Los U ^nlos U T^out- Since ^los^ ^nlos and 
T^out are independent, the density (or intensity), A( ), of L 
is equal to the summation of the intensities of ^los. ^nlos 
and ^ouT, i-e., A([0,x)) = Alos ([0, x)) + Anlos ([0:^;)) + 
Aout([0,x)). Since the path-loss of the links in outage is 
infinite, by definition Aout ([0?^)) = 0- On the other hand. 
As (•) for s = {LOS,NLOS} can be computed by using 
mathematical steps similar to the proof of [13, Lemma 1]. 
More specifically, we have; 

p-\-CO 

As{[0,x)) = 2 tt\ J n yx - {Ksrf‘fs{r)rdr (15) 

Equation (|9]l follows by inserting ps (•) of Q in (ITSI l and 
by computing the integrals with the aid of the notable result 

e-<=^rdr = (l/c^) (e”™ -f ae"™ - 6"°'' - 

B. Proof of Proposition Q] 

Erom © and the coverage probability can be formu¬ 
lated, by definition, as shown in (IThT l. Let us denote the first 
and second addends in ([T6] | by (•), where s = LOS and 

s = NLOS, respectively. Each addend can be computed by 
using the following results; 
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